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ABSTRACT 

The recent HEGRA observations of the blazar Mkn 501 show strong curvature in 
the very high energy 7-ray spectrum. Applying the 7-ray opacity derived from an 
empirically based model of the intergalactic infrared background radiation field (IIRF), 
to these observations, we find that the intrinsic spectrum of this source is consistent 
with a power-law: dN^/dE oc with a = 2.00 ± 0.03 over the range 500 GeV 
- 20 TeV. Within current synchrotron self-Compton scenarios, the fact that the TeV 
spectral energy distribution of Mkn 501 does not vary with luminosity, combined with 
the correlated, spectrally variable emission in X-rays, as observed by the BeppoSAX 
and RXTE instruments, also independently implies that the intrinsic spectrum must be 
close to a = 2. Thus, the observed curvature in the spectrum is most easily understood 
as resulting from intergalactic absorption. 

Subject headings: BL Lac objects: individual: Mkn 501: active galaxies - 7-rays: theory 
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1. Introduction 

Imaging air Cerenkov telescopes (lACTs) are 
currently able to detect 7-ray photons of TeV 
(10^^ eV) energy from BL Lac objects within one 
hour of observation and to measure their spectra 
using a few hours of a good data (for a review 
see Weekes et al. 1997 ). The objects of this class 
detected to date are, in order of increasing red- 
shift: Mkn 421 {z = 0.031) (Punch et al. [1992 



Petry et al. |1996|) , Mkn 501 {z = 0.034) (Quinn 
et al. [19961 Bradbury et al. [iM, lES 2344+514 
{z = 0.044) (Catanese et al. |1998|) and PKS 2155- 
304 {z = 0.117) (Chadwick et al. [l998[) . Measure- 



ments of the spectrum can be made over the en- 
ergy range from 200 GeV to about 10 TeV using 
this technique. Because these photons interact 
with infra-red radiation to form electron-positron 
pairs, the signal is expected to be attenuated by 
absorption both within the source itself and in 
the intergalactic medium. Thus, it is possible to 
use the observations to study the intergalactic 
infrared radiation field (IIRF), given some gen- 
eral, model-dependent constraints on the spec- 
trum intrinsic to the source (Stecker, De Jager, 
Salamon 1992[ ). Determining the IIRF, in turn, 
allows one to model the evolution of the galaxies 
which produce it. 

In this Letter we analyze recent observations 



of the object Mkn 501 (Konopelko et al. [1998b[) 
which are unique both for the quality of the 
spectra obtained and their energy range (up to 
20 TeV). These data, taken during a period in 
which the intensity of the source varied strongly, 
show a pronounced curvature in the spectrum, 
being significantly softer (steeper) towards higher 
energy. We unfold these data using the upper 
curve for the spectral energy distribution of the 
IIRF given by Malkan & Stecker ( [199^ , hence- 
forth MS98), with the corresponding 7-ray opac- 



ity as calculated by Stecker & De Jager ( 1998 , 
henceforth SD98) and find that the intrinsic spec- 
trum is flat, dN/dE oc E~'^. The implications of 
this result for both the absorption model and the 



synchro-self-Compton emission model are dis- 
cussed. 

2. Measurement of spectrum of Mkn 501 

The BL Lac object Mkn 501 exhibited strong 
emission in TeV 7-rays from March to October, 



1997 (Protheroe et al. [1998^ . During this pe- 
riod the source was continuously monitored by 
several ground-based imaging air Cerenkov tele- 
scopes, including the HEGRA stereoscopic sys- 
tem of 4 imaging air Cerenkov telescopes (Aha- 



ronian et al. 1997a), which observed it for a 
total exposure time of 110 hours (Aharonian et 



al. 1999). The unprecedented statistics of about 
38,000 TeV photons, combined with the good 
energy resolution of ~ 20% over the entire en- 
ergy range and with detailed studies of the de- 
tector performance (Konopelko et al. 1999| ), al- 
lowed a determination of the spectrum in the en- 
ergy range 500 GeV to 24 TeV (Konopelko et 
al. 1998b[ ). The Mkn 501 energy spectrum mea- 
sured by the HEGRA collaboration extends well 
beyond 10 TeV, where uncertainties related to 
the saturation effect could in principle play a role. 
Various data consistency checks were performed 
in order to avoid these effects. In addition, si- 
multaneous observations of the Crab Nebula (the 
standard-candle TeV source) were undertaken 
from 1997 September to 1998 March. Using sim- 
ilar data analysis, the Crab Nebula spectrum de- 
rived from the HEGRA data was found to be a 
pure power law with a differential spectrum index 
of 2.6 over the energy range 500 GeV - 23 TeV 
(Konopelko et al. 1998a[ ). These results are con- 
sistent with previous measurements by the Can- 
garoo group in the energy range 7~50 TeV (Tan- 



imori et al. 199< 



The flux of 7-rays, from Mkn 501, averaged 
over the entire observation period, was about 
three times that of the Crab Nebula (3 "Crab"). 
Averaged over each day, the 7-ray rate showed 
strong variations, with a maximum of 10 Crab 
detected on 26/27 June 1997. As remarked by 
Aharonian et al. ( [1997b| ), the hardness ratio of 
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the steepening Mkn 501 spectrum appears to be 
independent of the absolute flux. The high 7- 
ray detection rate provided event statistics of a 
few hundreds within 1 day's observations (~ 3— 5 
hours) which suffices to evaluate the energy spec- 
trum over the range 1-10 TeV. The analysis of 
the spectral shape on a daily basis did not re- 
veal any substantial correlation between the 7- 
ray flux and the spectral behavior (Aharonian 
et al. 1999| ). This justifies the presentation of 
a time-averaged energy spectrum of Mkn 501 in 
its active state, which is shown in Figure 1 over 
the energy range from 500 GeV to 24 TeV. The 
vertical error bars in this figure correspond to 
statistical errors. Note that the systematic er- 
rors at energies below 1 TeV appear to be quite 
large, reaching ~50 % at 500 GeV. 

Over the entire range, the spectrum shows a 
gradual softening towards higher energy. The 
19-24 TeV energy bin contains a signal with a 
significance of 3.7a. However, the steep energy 
spectrum and 20 % energy resolution do not per- 
mit one to exclude the interpretation that these 
7-rays may have spilled over from the lower en- 
ergy bins. Therefore, the energy spectrum is con- 
sistent with the hypothesis of a maximum energy 
for the detected 7-rays of ~ 18 TeV. The shape of 
the energy spectrum is well described by a power- 
law with an exponential cutoff. A fit of the data 
over the energy region where the systematic er- 
rors are small, i.e., from 1 TeV to 24 TeV, gives 



dN/dE 
A 
a 
Eo 



AE^'' exp i-E/Eo) [cm'^s-^TeV"^] . 
(9.7 ± 0.3 (stat) ± 2.0 (syst)) • 10"" 
-1.9 ±0.05 (stat) ±0.05 (syst) 
5.7 ± l.l(stat) ± 0.6 (syst)TeV 

(1) 



The logarithmic slope of the energy spectrum 
("power-law index") is 1.8 in the energy range 
1-5 TeV, and 3.7 above 5 TeV. 

Independent measurements of the Mkn 501 
TeV energy spectrum by the Whipple Observa- 



agreement with these results. The best fit to the 
data presented by the Whipple group agrees pre- 
cisely with a fit to the HEGRA data in the energy 
range 500 GeV - 10 TeV. However the HEGRA 
group measured the spectrum well above 10 TeV 
where it exhibits a further steepening. 

3. Absorption on the diffuse intergalactic 
infra-red background 

The formulae relevant to absorption calcula- 
tions involving pair-production are given and dis- 
cussed in Stecker, De Jager & Salamon ( |l992| ). 
For 7-rays in the TeV energy range interacting 
at redshifts z <C 1, the pair-production cross sec- 
tion is maximized when the soft photon energy 
is in the infra-red range: 



X{E^) ~ A, 



E.y 



2.4E^^TeV /xm (2) 



tory (Samuelson et al. 1998| ) are in very good 



where Ae = h/{mec) is the Compton wavelength 
of the electron. For a 10 TeV 7-ray this cor- 
responds to a soft photon in the mid infra-red 
region of the spectrum, having a wavelength 
around 24 /um. Pair-production interactions take 
place with photons over a range of wavelengths 
around the optimal value, as determined by the 
energy dependence of the cross section. Stecker 
& De Jager (SD98) have computed the absorp- 
tion coefficient of intergalactic space using a new, 
empirically based calculation of the spectral en- 
ergy distribution (SED) of intergalactic low en- 
ergy photons (MS98). Assuming that the IIRF 
is basically in place by a redshift ~ 0.3, hav- 
ing been produced primarily at higher redshifts 
(Stecker & De Jager |l99^ ; Madau |T995| ), 

SD98 limited their calculations to z < 0.3. Evo- 
lution in stellar emissivity affects the predicted 
IIRF and is expected to level off or decrease at 
redshifts greater than ~ 1.5 (Madau |1996| ). In 
this paper, we assume that evolution continues 
up to z = 2, leading to the higher of the two 
IIRF used by SD98. This is more consistent 
with recent data on IR galaxy evolution, dust 
absorption, and the lower limits from IR galaxy 
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counts (Stecker |1999| ). To compute the absorp- 
tion, we adopt the SD98 parametric expressions 
for t{E,z) for z < 0.3, taking a Hubble constant 
of Ho = 65 km s'^Mpc'^ 

The unfolded HEGRA data points are also 
shown in Figure 1, together with a fit to a power 
law energy spectrum. We find 



dN^/dE 



1.32 ±0.04 • 10"^° X 
(i?/lTeV)-2-00±0-03 

photons cm^^ s~^ (3) 



with = 18.6 for 15 degrees of freedom, giv- 
ing a high chance probability of 0.2. In the 
mid-energy region 1-10 TeV, where the measured 
energy spectrum is very well-defined, the data 
points deviate from the fit by less then 15%, 
which equals the estimated systematic error. Note 
that both the statistical and the systematic un- 
certainties of the spectrum measurements in- 
crease towards the upper end of the energy range, 
where they reach 30% and 60%, respectively. 
Thus, the data points of the unfolded spectrum 
are consistent, within the statistical and system- 
atic errors, with the simple power law fit of dif- 
ferential spectrum index 2.0. Our analysis shows 
that fitting the unfolded spectrum using an ad- 
ditional exponential term, result in a x^-value of 
the same magnitude. 

4. Implications for the intrinsic spectrum 
and the inter-galactic absorption 

In order to understand the implications of 
the absorbed and de-absorbed spectra shown in 
Figure 1, it is necessary to adopt a model or 
scenario for the production of TeV photons in 
the source. Of the many suggestions in the lit- 
erature, interest has recently centred on those 
in which a single population of relativistic elec- 
trons is responsible for both the TeV photons 
and for photons in the X-ray region of the spec- 
trum, as in the synchrotron self-Compton (SSC) 
and external Compton models (e.g., Bloom & 
Marscher fmE[ , Inoue & Takahara |l99|; Ghis- 
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Fig. 1. — The energy spectrum of Mkn 501 as 
measured by the HEGRA lACT array (open circles) 
(Konopelko et al. 1998b). The combined power law 
plus exponent fit of the HEGRA data is shown by the 
dotted-dashed curve. The Mkn 501 spectrum mea- 
sured by the Whipple group (filled circles) is taken 
from Samuelson et al. ( |1998| ). Also shown are the de- 
absorped HEGRA points (open circles) found using 
the optical d epths to absorption calculated by Stecker 
& De Jager ( 1998 ) for the "high" intergalactic infra- 
red radiation field. The power-law fit to these data 
is shown by the solid line, and has a photon index of 
a = 2.00 ±0.03. 



ellini & Madau 1996 ; Dermer, Sturner & Schlick- 
eiser p97| ; Mastichiadis & Kirk p97| (henceforth 
MK97); Sikora et al. 1997 ; Georganopoulos & 



Marscher 1998; Ghisellini et al. 1998; Levinson 



1998). These models are favored because they 
provide a natural way to understand the simi- 
lar variability timescales of the X-ray and TeV 



emission. 



During the period of the TeV observations, 
Mkn 501 was also observed in X-rays using the 
BeppoSAX instrument (Plan et al. 1997] ) and the 
Rossi X-ray Timing Explorer (Lamer & Wag- 
The spectrum in this energy range 



ner 



19981) . 



varied strongly, with generally a very hard spec- 
tral index extending to much higher energies 
^ 100 keV than during less active epochs. There 
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was a close temporal correlation between the X- 
ray and TeV fluxes, further strengthening the 
case for the origin of X-rays and TeV in a com- 
mon electron population. 

Extensive studies of the variability properties 
of the SSC model have been undertaken (MK97), 
which show that the mechanism most likely to be 
responsible for the variability shown by the ob- 
ject Mkn 421 is a change in the maximum energy 
(expressed Lorentz factor: 7max) to which 
electrons are accelerated. These results were also 
apphed to Mkn 501 (Mastichiadis k Kirk |1999| , 
henceforth MK99), where, once again, variabil- 
ity induced by a change in 7max appears to give 
a reasonable fit to the X-rays and to the TeV 
data then available. 

For such models, the most important property 
revealed by the HEGRA observations discussed 
above is the lack of variation in the TeV spec- 
trum, despite the fact that the correlated vari- 
ations in the X-rays show strong spectral vari- 
ations, consistent with an increase in 7max- In- 
spection of the model light curves in the TeV 
range show that the softer the spectral slope be- 
comes, the more sensitively it reacts to changes 
in 7max- This is because intrinsic spectra softer 
than a photon index of a = 3 are a direct result of 
the electron cut-off, whereas harder spectra (pho- 
ton index a ~ 2) can be formed by power-law 
electrons scattering off a range of target photon 
energies. 

Thus, on the basis of the homogeneous SSC 
model, we find that intensity variations with con- 
stant spectral shape imply an intrinsic spectrum 
with a photon index of a ~ 2. This provides ad- 
ditional evidence that the spectrum of Mkn 501, 
which has a slope a ^ 3.7 above 5 TeV, must be 
modified by the effects of inter-galactic absorp- 
tion. Of the two IIRF considered by SD98, only 
the higher provides sufficient absorption to ac- 
count for such a strong modification. We note 
again that lower limits from galaxy counts in 
the mid-IR, as well as other observational data 
(Dwek et al. |199|), favor the high IIRF (Stecker 



19991) . 

As well as constraining the intergalactic ab- 
sorption, the observation of an intrinsic spectrum 
of a ~ 2 at 10 TeV requires a higher Doppler fac- 
tor than considered by MK97 and MK99. Using 
the approximate scaling laws presented in MK97, 
we can estimate that a Doppler boosting factor 
of ~ 50 suffices to produce a 2 at 10 TeV, 
and have confirmed this by running full simula- 
tions. It is interesting to note that such boost- 
ing factors, although larger than values measured 
in sources which display apparent superluminal 



motion (Vermeulen & Cohen 1994), seem to be 
indicated both by observations of intra-day vari- 
ability (Wagner & Witzel 1995| ) and of extremely 
rapid variations in the TeV flux of blazars (e.g., 
Gaidos et al. [l996| ). 

A completely model independent conclusion is 
still elusive, and will remain so until observations 
of comparable quality on other blazars of differ- 
ent redshifts are available. Nevertheless, the two 
independent arguments we have presented favor- 
ing an intrinsic emission spectrum close to a = 2 
indicate that the effect of absorption by the in- 
tergalactic infra-red background radiation in the 
spectrum of Mkn 501 is strong, and suggest a 
Doppler boosting factor for this source of ^ 50. 
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